arXiv:astro-ph/9701136vl 20 Jan 1997 


OU-TAP-53 
January 1997 


STABILITY OF SCALE-INVARIANT COSMOLOGICAL 
CORRELATION FUNCTIONS IN THE STRONGLY 
NON-LINEAR CLUSTERING REGIME 

Taihei Yano and Naoteru Gouda 

Department of Earth and Space Science, Graduate School of Science, Osaka University 

Toyonaka, Osaka 560, Japan 
E-mail: yano, gouda@vega.ess.sci.osaka-u.ac.jp 

ABSTRACT 

We investigate stability of the scale-invariant solutions of the BBGKY 
equations for two-point spatial correlation functions of the density fluctuations 
in the strongly non-linear regime. In the case that the background skewness of 
the velocity held is equal to 0, we found that there is no local instability in the 
strongly non-linear regime. The perturbation does not grow nor does it decay. 

It has an only marginal stable mode. This result means that no special value 
of the power index of the two-point spatial correlation function are favored in 
terms of the stability of the solutions. In other words, the argument about the 
stability does not determine the power index of the two-point spatial correlation 
functions in the strongly non-linear regime. 

Subject headings: cosmology:theory-large scale structures-correlation function 


1. INTRODUCTION 

The large scale structure formation is one of the most important problems in 
the cosmology. It is generally believed that these structures have been formed due to 
gravitational instability. Hence it is very important to clarify the evolution of density 
huctuations by gravitational instability. Here we consider the density huctuations of the 
collisionless particles such as dark matters because our interest is concentrated on the effect 
of the self-gravity. 
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When the density fluctuations are much smaller than unity, that is, in the linear 
regime, time evolutions of the small density fluctuations can be analyzed by making use of 
the linear theory. In this regime, we can understand analytically how the small fluctuations 
grow(Peebles 1980, 1993). But when the density fluctuations are much larger than unity, 
that is, in the strongly nonlinear regime, the analytical approach is very difficult. But we 
believe that it is very necessary to understand clearly the nonlinear behavior of the density 
fluctuations on these strongly nonlinear scales because galaxy formations are much related 
to the density fluctuation on these small scales and also it is a very interesting academic 
problem for the nonlinear dynamics of the self-gravity. 

As one way to quantify the clustering pattern, we generally use a two-point spatial 
correlation function of density fluctuations. In the strongly non-linear regime, it is found 
from iV-body simulations that two-point spatial correlation functions obey the power law. 
This result is reasonable because the self-gravity is scale-free. Two-point spatial correlation 
function has been investigated by various methods because the power index of the two-point 
spatial correlation function is a good indicator representing the nonlinear dynamics of the 
self-gravity in this regime. 

The power index of the two-point spatial correlation function have been usually 
analyzed by iV-body simulations(Frenk, White, & Davis 1983; Davis et.al. 1985; Suto 1993 
and references therein). However, the physical process that determine the value of the 
power index cannot be clarrifled only by the iV-body simulations. There are some other 
methods besides this numerical simulation. One of them is the analysis by the BBGKY 
equations. The work by Davis & Peebles (1977; hereafter DP) is a pioneer for the analysis 
by the BBGKY equations. They showed the existence of the self-similar solutions for spatial 
correlation functions under some assumptions. Then it is shown that the power index 7 of 
the two-point spatial correlation function in the strongly nonlinear regime is related to the 
initial power index n of the initial power spectrum P{k) as follows; 

^(r) oc r-^(^ > 1 : 7 = ^^^^) ( 1 ) 

One of the assumptions that DP adopted is called the stable condition. This condition 
is that the mean relative physical velocity in the strongly nonlinear regime is equal to zero. 
This condition was tested by the iV-body simulations (Efstathiou et al. 1988; Jain 1995), 
but this condition is not completely verified. Furthermore the stability of the scale-invariant 
solutions in the strongly nonlinear regime, which DP derived, was investigated by the linear 
perturbation theory(Ruamsuwan & Fry 1992;hereafter RF) because there is no guarantee 
that such solutions can exist stably. And it is found that the solutions are marginally stable. 
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As for the physical process determining the power index in the strongly nonlinear 
regime, there are other analysis besides one proposed by DP. One of them is given by 
Saslaw(1980). He concluded that the power index 7 approaches to 2 by using the cosmic 
energy equation under some assumptions while some numerical simulations do not support 
this result (Frenk, White & Davis 1983; Davis et ah 1985; Fry & Melott 1985). 

There is another idea as follows; when the initial power spectrum has the sharp cut-off 
or the initial power spectrum is scale-free with negative and small initial power index, 
then there appear anywhere caustics of the density fields. In these cases, the power index 
is irrespective of the detailed initial conditions after the first appearance of caustics on 
the small scales around the typical size of the thickness of caustics(in three-dimensional 
systems, they correspond to the pancake structures of highly clustered matters). The power 
index is determined by the type of these caustics which is classihed in accord with the 
catastrophe theory. This idea is verihed in the one-dimensional system (Kotok & Shandarin 
1988; Gouda & Nakamura 1988, 1989), the spherically symmetric systems(Gouda 1989), 
the two-dimensional systems (Gouda 1996a) and also the three-dimensional systems(Gouda 
1996b). In these cases, it is suggested that 7 0 on the small scales. As we can see 

from the above arguments, we believe that there are still uncertainties about the physical 
processes which determine the value of the power index. 

Yano & Gouda (1996;hereafter YG) investigated the conditions that determine the 
power index of the two-point spatial correlation function in the strongly nonlinear regime 
by analysing the scale-invariant solutions of the BBGKY equations in this regime. YG does 
not adopt the following assumptions that DP adopted; (1) the skewness is equal to 0, (2) 
three-point spatial correlation function is represented by a product of the two-point spatial 
correlation function, (3) the stable condition is satished. As a result, YG obtained the 
relation between a mean relative pecurier velocity, skewness, three-body correlation function 
and the power index of the two-point spatial correlation function. YG found that the stable 
condition is not independent one of the other assumptions (1) and (2). That is, the stable 
condition is satisfied only when both assumptions (1) and (2) are satished. The assumptions 
( 1 ) and ( 2 ) cannot be generally satished in all cases and then there is no guarantee that 
the stable condition are correct. Furthermore YG suggests from the physical point of view 
the probable range of the mean peculier velocity which includes the the value given by the 
stable condition. This fact results in the possibility that the power index of the two-point 
spatial correlation functions takes various values according to the mean peculier velocity. 
Indeed, YG showed that the mean relative physical pecurier velocity can take the value 
between 0 (stable clustering) and the Hubble expansion velocity (comoving clustering). 
When the stable clustering picture and the self-similerity are satisfied, the power index of 
the two-point spatial correlation function have the value that DP derived. On the other 
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hand, when the comoving clustering picture and the self-similerity are satished, the power 
index of the two-point spatial correlation function have the value of 0. And this value is 
consistent with the result from the catastrophy theory(Gouda& Nakamura 1988 1989;Gouda 
1989 1996a;b). Although we found that there exist various scale-invariant solutions with 
different power index, whether these solutions are stable or not is an another interesting 
probrem. Some values of the power index of the two point spatial correlation function may 
not be able to be taken in the real world if these solutions are unstable. RF investigated the 
stability of the DP solutions by making use of the linear perturbation. And they showed 
that the perturbations of the solutions is marginal stable. But they investigated only in the 
DP case and did not investigate the other solutions that YG obtained. Furthermore, as we 
will discuss later, RF mistook the way of providing the perturbation about the skewness 
although their result is correct fortunately due to the reason which we will show later 
(§3). Furthermore they did not comment about the ’strange’ growing mode which exist in 
their solutions of the perturbation. We will show that this ’strange’ modes are resulted 
from the fact that RF put the inapplicable form of the perturbation which diverges on 
small or large scales. So, in this paper, we investigate the stability of the general solutions 
that YG obtained by putting the applicable form of the perturbation. We will derive the 
perturbation equation by perturbed BBGKY equation from the background scale-invariant 
solutions in the strongly nonlinear regime. And we will show the solution of the linear 
perturbation equation when we assume the appropriate form of the perturbation which is 
well-dehned. 

In §2, we briefly show the BBGKY equations that we use in this paper and also show 
the scale-invariant solutions in the strongly non-linear regime that YG obtained. In §3, we 
consider the stability of these solutions of the BBGKY equations by the analysis of the 
linear perturbation. Finally, we devote §4 to conclusions and discussions. 


2. BASIC EQUATIONS AND THE SCALE-INVARIANT SOLUTIONS 

At hrst, we show the BBGKY equations that we use in this paper. These equations 
are time evolution equations of the statistical value such as two-body correlation function, 
three-body correlation function, and so on. These equations can be derived from the 
ensumble mean of the Vlasov equation (DP,RF). The A-th BBGKY equation represents 
the time evolution of the A-body correlation function. We are now interested in the 
two-body correlation function, and then, we use the second BBGKY equation by taking 
the momentum moment of this equation. Zeroth moment of the second BBGKY equation 
becomes the time evolution equation of the two-point spatial correlation function. This 
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equation involves the first moment term, that is, the term of the mean relative peculier 
velocity. The hrst moment of the second BBGKY equation becomes the time evolution 
equation of the mean relative peculier velocity. This equation involves the second moment 
term of the relative peculier velocity dispersion If, E, where. If ,and E are pararel 
and transverse component of the mean relative peculier velocity dispersion, respectively 
(DP,RF,YG). The second moment of the second BBGKY equation becomes the time 
evolution equations of the relative peculier velocity dispersion. These equations invole the 
skewness of the velocity held in the same way. In deriving the BBGKY equations, DP 
assumed that the skewness is equal to 0, and also the stable condition in which the mean 
relative peculier physical velocity is equal to 0 ((u) = —ax). Furthermore DP assumed that 
a three-point spatial correlation function can be represented by a product of the two-point 
spatial correlation function as follows; 


Cl23 — <5 ('Ci2'C23 + ^ 23^31 + ^ 31 ^ 12 )- (2) 

On the other hand, RF incorporated the skewness in the equations. RF expressed the 
skewness by using the values A, and B. These are related to our form by the next relations. 


R(3unRi7’-t-= 3(u)n-|-S||, RuE = (u)E-f s_L, (3) 

where v in RF is the same as (v) in our paper, and E in RF is also the same as our E. The 
dehnition of IlijiT’ is different from ours (IlRiT’ -(- = If). Here sy, s_l are, as we will dehne 

later, the parallel component and the transverse component of the skewness, respectively. 
When we investigate the scale-invariant solutions of the BBGKY equations, the difference 
of the dehnition of the valuables in RF form and ours is not important. But when we 
perturbe the each valuable, the skewness also must be perturbed independently. However 
it must be noted that RF treated A and B as constant values. Then RF did not correctly 
treat with the perturbation of the skewness. RF used the same assumption that DP used 
about the three-point spatial correlation function. Furthermore, they also used the stable 
condition. We do not know whether the assumption of the three-point spatial correlation 
function is correct or not. Then we assume the three-point spatial correlation function by 
the following form; 


Cl23 = Q(ei'2+'e23^' + e23^'e3i+' + , (4) 

where 5 is a constant value. This is not a general form, but an extention of the form that 
DP adopted. We use this form as a preliminary step for our analysis. Furthermore the 
stable condition is a special case of the mean relative peculier velocity. Hence we do not 
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assume the stable condition and consider the general condition as shown in YG. We are 
interested in the strongly nonlinear regime, and then we take the non-linear approximation. 
In the strongly nonlinear regime ( x -C 1), the two-point spatial correlation function is 
much larger than unity, 3> 1. In this limit, we obtain the following four equations (DP 
1977;RF 1992;YG 1996); 

11 ^ 

-—-I-—= 0, (0th moment) (5) 

ot a ox 

ax 

2GmnaQ— [ + ^(^)^^^}^(^ - x^+^d^Xs = 0, 

'T* / ^ 

Jb J *^31 

(1st moment) (6) 


ax^ dx 


>"en) - 


dt 


1 d r 


f>l3 rf'y 


4 e 


+ AGmnaQ*^^-^ f + C(2:)^’^‘^}^(2: - x)^+‘^(n'^)d^X3 = 0, 

X J ^ 3 ]^ 

(2nd moment : contraction 1) 




+ -^7?-+ s_l }1 = 0 , 


ax* dx 1 


(2nd moment : contraction 2) 


(7) 


( 8 ) 


where G is the gravitational constant, m is the mass of a particle, n is the mean number 
density of the particles, Q* is the coehcient of the hrst momentum moment of the three-body 
correlation function (RF,YG). As YG commented, the fourth term of eq.(|^ is not satished 
in general although this term is correct in the strongly non-linear regime in general. 
However DP showed the existence of the three-body correlation function which gives this 
term. As commented later, the condition that Q* = Q is consistent with the zero skewness 
case [see eq.(^^]. Since we consider the case of the zero skewness, we assume the fourth 
term of eq.(|^) is satished in our analysis. 


We express the skewness by the following expression. 


= ((!. - (!.))“(.« - (.«»'’(.« - = S|| + Si Pplf, 


pafS'y 

PPP 


x^ x^ x'^ 


x^ 


rviOL rfPf 0 rt^'y 


X 


X 


X 


X-" 


(9) 


( 10 ) 
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where the subscripts p and t represent the parallel and transverse component of each 


two particles, respectively. Here Pppt and Put vanish because of the symmetry of the 
background universe. We use the Einstein-de Sitter Universe through this paper because 
we are interested only in the scale-invariant correlations. 

Here we consider the scale-invariant solutions of these equations. In the strongly 
nonlinear regime, it is naturally expected that the effect of the nonlinear gravitational 
clustering dominates and then the solutions in this regime have no characteristic scales, 
that is, they are expected to obey the power law due to the scale-free of the gravity. Then 
we investigate the power law solutions of the (u),n and S (YG). We assume that the 
two-point spatial correlation function ^ is given by 




( 11 ) 


Then, we obtain from the dimensional analysis in eq.(^ 


{v) = —hdx, 

= (3 - 7)/i. 


( 12 ) 


In this case, the solutions of the other valuables are given by 


s 


H 


Eoa^('+25)-1^2-^(l+25)^ 

S||„aa^(i+25)-i^3-7(i+2^)^ 


X' 


(13) 


From eq. (j^) , it is found that 


2/3(1 -|- 5) -[-1 — {7 — 2^(1 -|- 5)}(h — A) — 0, 


(14) 


and then 


1 + {7-27(1 + 5)}A 


(15) 


1-6(5 



( 16 ) 



As we can see from eq.(^), the parameter h can take varions valnes according to the 
skewness A, the power index 7 and 5. Only when A = <5 = 0 is satished, the stable condition 
{h = 1) is correct. If the similarity solntions exist, the power index of the two-point spatial 
correlation fnnction can be represented by the following form (Padmanabhan 1995;YG 
1996); 


3(3 -|- n)h 
2 -|- (3 -|- n)h 


( 17 ) 


Then, even when we assnme that the self-similarity solntions exist, the power index of the 
two-point spatial correltion fnnction can take varions valnes according to the mean relative 
pecnlier velocity, h (if the stable condition is satished, i.e., h = 1, the resnlt of DP[eq.(2)] is 
reprodnced). 


3. LINEAR STABILITY OF SCALE-INVARIANT SOLUTION 

In the previons section, we showed that there are varions scale-invariant solntions in 
addition to the solntions that DP obtained (YG 1996). Bnt there is no gnarantee that all 
the solntions can exist stably. As RF did, we investigate the stability of the solntions in 
the strongly nonlinear regime by making nse of the linear pertnrbations theory. Now we 
pertnrbe the two-point spatial correlation fnnction from the scale-invariant solntion that we 
obtained in the previons section as follows; 

e' = e(l + A^), ( 18 ) 

where ^ is the scale-invariant solntion, is the pertnrbed one and A^ -C 1. We also 
pertnrbe the other valnables snch as the mean relative pecnrier velocity (n), the relative 
pecnrier velocity dispersions 11 , S and the skewness sy, sj_ in the same way. 

The eqnations of the linear pertnrbations are following; 

^ moment) (19) 

Y||,.^n{A, + An}| - + 

+ 2GmnaQa;f (1 + = 0, 

( 1 st moment) 


( 20 ) 
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l|[a^5nl{A, + An} 


ax^ dx - 

4^ 


+ An} + S||{Aj + A^,}} 


-{(v)S{A^ + A(„) + As} + s±{A^ + A^j^}} 

CLOC 

+ AGmnaQ*x{v)(^^^^^\{l + + M^(^i+s)^{v)] = 0 , 

(2nd moment : contraction 1) 


( 21 ) 


j__a 

dt 

+ 


[o‘" 5E{A5 + As}] 


1 d r 

ax^ dx - 


+ ‘^{v) + As} + s_l{A^ + As^}] 


(2nd moment : contraction 2) 


0 . 


( 22 ) 


Here it should be noted that in deriving the above equations (|1^)-(|^) we neglect some 
terms in full equations because of the following two reasons; one of them is that some terms 
are high order ones in the strongly non-linear limit. Another reason is that some terms 
are the higher order ones in the limit of small perturbation (larger than the hrst order 
perturbation), he ordering parameters in bothlimits are generally independent of each 
other. However we consider the case that the higher order terms in the strongly non-linear 
limit are much smaller than the hrst order terms of the perturbation. Then we can derive 
the above equations by neglecting higher order terms in the nonlinear approximation. This 
case means that we neglect the effects of the higher order terms in the strongly non-linear 
approximation in considering the linear perturbation. 


RF expected the power law perturbation given by the following; 

Ag = e^aFx^^. 


(23) 


In this case, and ^ are given by 


^7(1+<5) 


/i 


J 

J 


(y (y 

X Z Xsi 

X Z ^ ^ X2l 


- s - — 
x’ X21 


(i + y-29^)‘''7 


(24) 

(25) 

(26) 


The integrals should not diverge. Then, 2 — 7(1 -f- 5) > 0 must be satished for ?/ —> 0 
and 7(1 -|- 5) > 0 for 1/ —>• cx) in the M. Furthermore, 2 -f- g — 7(1 -|- 5) > 0 must be satished 
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for 1 / —>■ 0 and g — 7(1 + 5) <0 for cxo in the M'. As a resnlt, the following relations 
must be satished; 

0<7(1 + 5)<2, 7(1 + 5)-2 < g < 7(1 + 5). (27) 

The four perturbation equations ( 0 )-(S) are a little different from the equations that 
RF derived. This is because in perturbating the three-body correlation term, RF devide 
the {v 2 i) into (^ 23 ) + ('t'si) artihcially and perturbed {V 23 ) and (^ 31 ) independently, which 
is incorrect treatment. Furthermore, although they obtained the g-dependent of M', they 
used the value of M' only at g = 0 in solving the perturbation equations. 

Here we consider the following form of the perturbation. We need to put the 
perturbations that do not diverge in the strongly non-linear limit {x —> 0 ) and the linear 
limit {x 00 ). Then we put the perturbation of the two-point spatial correlation function 
as follows; 


= e^a^x‘', (28) 

g = |g|i, (29) 

where |g| is the real number, that is, g is the pure imaginary number while RF adopted 
the real number. In this case, the perturbations do not diverge in the strongly non-linear 
limit and also in the linear limit. If g is real number as RF adopted, the perturbation 
diverges on some scales. When g is negative, the perturbation diverges in the non-linear 
limit (a: —> 0). On the other hand, when g is positive, the perturbation diverges in the linear 
limit {x —> cxo). The perturbations of the other valuables also can be put in the same form. 
From the dimensional analysis, all perturbations must have the same value of the power 
(g and p) as that of the two-point spatial correlation function. When g is pure imaginary 
number, Is”^! = 1 and = 1 . And the following relations are satished; 


l^7(l+<5),ql < M^(l+5),q=0 — 2My(l+5). 


(30) 


When 0 < 7(1 + 5) < 2, both M^(^i+s) and are hnite for any value of pure imaginary 

number g. 


These perturbed BBGKY equations are not closed by themselves in general and higher 
moment equations are needed. But if the coefficient of the perturbation of the skewness 
(A^ii and A^^) are equal to 0 , the perturbation equations for the other valuables have no 
relation to the perturbation of the skewness. In this case we can solve these perturbation 
equations independently from the higher moment equations. As we can see from eqs.(^) 
and (^), when the skewness of the background velocity held is equal to 0 , the coefficient 
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of the perturbation of the skewness becomes 0. Here we consider only this case that the 
skewness is equal to 0. Then we can neglect the higher moment equations. RF treated A 
and R as a constant value. So, the treatment of the perturbation of the skewness itself 
was mistaken. But when the skewness is equal to 0, RF treatment happen to be correct 
fortunately. In this case, the above equations are rewriten by putting the perturbations 
which have the power law given by eqs. 


and (^9]) as follows; 


(p - hq)A^ - h(3 - 7 + g) = 0, 


(31) 


[2{7(1 + 5) — 2 + cr}(l + 2^) + g + 2(1 + 6)Dkq\A^ 

+ {4 - 27(1 + 6 ) + g}An - 2aAs = 0 , (32) 


[1 — 15h + (6 + 47)h(l + 5) + 4ha + p — 3hq 

-8h^{7(l + 5) - 2 + a}(l + 5) - 4h^(l + 6)Dkq]A^ 

+ [—3h{5 — 27(1 + 5) + g} + 4ha — 4h-—{7(1 + 5) — 2 + a}] A^^,^ 

V 

+ [1 — 15h + (6 + 47)h(l + 5) + p — 3hg] An + 4/iaAs = 0, (33) 


(1 — h + 6h6 + p — hq)A^ — h{7 — 27(1 + 5) + g} A(^) + (1 — h + 6h6 + p — hg)As = 0, (34) 


where 


= 7(1 + 5) - 2 + a, 

It is difficult to treat exactly M' as a function of g. We approximate M' by linear 
function of g. When g = 0, M'/M is equal to 2. Then we use the following approximation 
of M'/M in the above equations (|^) and (PI); 


M' 

H 


= 2 + fcg. 


(36) 
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where 

/ logs + \ogy]d^y 

The integral + y~'^^^~^^'>)(Py is dominated arround the s and y ^ 1. 

Arround the s and ?/ ~ 1, logs and \ogy have the values of order 1 . So, k has the value of 
order 1 . 

Furthermore we use the following relation that can be derived from the hrst moment 
equation [eq.(|D]; 


cl+25 

4 — 27(1 + (5) — 2(7 + 2Gmna? Qx^ = 0. (38) 

These four equations [eqs.(pTD-(p4D] can be rewriten by using a matrix expression. 


= (^5) ^n, As)- 


(39) 

(40) 


If there exist a non trivial solution, the determinant of the matrix Nij should be equal 
to 0. Now we consider the zero skewness case, there is the relation between the three-point 
spatial correlation function and the mean relative pecurier velocity by using eq. (P!5|) . 


5 


h-1 

6h 


(41) 


By using this relation, we can eliminate 6 in eqs. (|3^) - (^^ . Furthermore, from the hrst 
moment equation(^ and the second moment (contraction 1 ) equation (|^, we obtain the 
following relation; 


Cm Ac I 

(-h + l + 6h6)U + {5 - 27(1 + 5)}^-^ - 4Gmna^x^M^n+s)h(Q* - Q) = 0. (42) 

ax ax 

As we can see, when the skewness is equal to 0, Q* — Q must be 0 , that is, Q*/Q is equal 
to 1 . 


Then the determinant of N is given by 

det N =^(p- hqffiq, 7'. h, a, kD), 


(43) 


where 
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f{q,'j',h,a,kD) = {9h + kD{7h — l)}q‘^ 

+ (4 - 27 ' + 29h - lOy/i -2a + 8ha)q 
+ kD{—3 + 21h — 6'y'h)q + 12 — 67 ' — 6a, (44) 


and 

7 '= 7 ( 1 + 5). (45) 

Here / is the quadratic equation of q. Now we investigate whether the equation / = 0 has 
real solutions or not. Since we do not know the value 7', h, a, and kD in the strongly 
non-linear regime, we treat these values as parameters. Here we consider the allowed 
range of these parameter. As we can see from eq. ([27|) , the parameter 7' must be satished 
0 < 7 ' < 2. Since we have treated as ^ ^ in the strongly non-linear regime, 

should be higher order than y This results in that 2(1 -|- 5) > 1 must be satished. In this 
case, h > j must be satished from eq. (i3)- YG showed the probable range of the mean 
relative peculier velocity and obtained that the mean relative physical pecurier velocity 
must have the value between 0 and the Hubble expansion velocity. This means that the 
parameter h(relative velocity parameter) have the value between 0 and 1. So in this case, 
the parameter h should be in the range | < h < 1. We do not know the range about the 
parameters a and kD. But the parameter a sould have an order of 1. Hence we investigate 
a in the range | < cr < 2. The parameter D and k also take an order of 1 as seen from the 
eqs. (|35[) and (|37|), respectively. So we investigate kD in the range — 1 < kD < 1. 


In the above probable value of the parameters, we can easily ascertain that / = 0 has 
real solutions. In other words, / = 0 does not have the solutions of imaginaly number. 
Since we consider the case that q is imaginaly number, / = 0 can not be satished. p — hq 
should be 0 to satisfy that the determinant of the matrix N is equal to 0. In this case, p 
has the following value 

p = hq = h\q\i. (46) 

This means that the perturbations do not grow. And the solutions are stable. 


Furthermore we consider the strict condition which determines the stability of the 
two-point spatial correlation function. In the strongly non-linear regime, the mean relative 
peculier velocity has the value, (n) = —hax according to the process of clustering. In this 
case, the scale of a^x for the two particles whose mean comoving distance is x does not 
change as we can see from the following relation; 


d 

dt 



a 


h-l 


{ax + hdx) 
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= ^{{v) + hax) 

= 0 (47) 

Then we should determine the stability of the two-point spatial correlation function at the 
fixed scale of a^x. The solutions of the perturbation are rewritten by 

= e^{a^xY 


This perturbation does not grow nor does it decay. At the fixed scale of a^x, the 
perturbation never even oscillate. This means that the perturbation is marginal stable. 
RF used the real number q in investigating the behaviour of the perturbation. As we can 
see from eg. (|48|) , the perturbation in the p — hq = 0 mode works well even when q is real 
number. That is, the perturbation is marginal stable in this mode. However there also exist 
a ’strange’ growing mode for the real number of q because / = 0 is satished in this case. 
RF did not comment sufficiently about this ’strange’ growing mode. 


When q is negative, the perturbation diverges in the non-linear limit (x —> 0). On the 
other hand, when q is positive, the perturbation diverges in the linear limit (x —> cx)). Then 
this perturbation is not adequate in investigating the local stability of the two-point spatial 
correlation function in the strongly nonlinear regime. 


4. RESULTS AND DISCUSSION 

In this paper, we investigated the stability of the scale-invariant solutions of the 
cosmological BBGKY equations in the strongly nonlinear regime in the case that the 
skewness of the velocity field is equal to zero. The reason why we consider the only case that 
the skewness vanishes is that perturbed BBGKY equations for (v), H and S can be closed 
independently of the higher moment perturbations. When the power law perturbations are 
put in the solutions as RF put, that is, when q is real number, the perturbation in the 
linear limit or the non-linear limit diverges. When q is negative, the perturbations diverge 
and do not work in the non-linear limit. So q should be positive. On the other hand, when 
q is positive, the perturbations work well in the non-linear regime. In the linear regime, 
however, those perturbations may diverge if the power law form of the pertubations are 
retained and so the form of the perturbations should be changed in order not to diverge 
in the linear limit. At this case, it is insufficient to solve the non-linear approximated 
equations because the informations about the evolutions on all scales are needed. Then 
we investigated only the local stability of the non-linear regime. In investigating the local 
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stability of the two-point spatial correlation function in the strongly non-linear regime, we 
should put the perturbation whose value is zero or smaller value on the scales except for the 
strongly non-linear scales than that on the strongly non-linear scalesthe investigating scale. 
That is, we should put the wave packet-like perturbation. In order to put such a wave 
packet-like perturbation, q must be imaginaly number. In this case, we found that there is 
no unstable mode in the strongly non-linear regime. It seems stable for any value of the 
power index of the two-point spatial correlation function in the strongly non-linear regime. 
However we do not know whether a grobal unstability exist or not because we consider only 
the local stability. It is certain that there is no local unstability in the strongly non-linear 
regime. So, in the strongly non-linear regime, the solutions is marginally stable, and it does 
not seem that the power index of the two-point spatial correlation function approaches to 
some stable point values. The power index of the two-point spatial correlation function 
that was derived by DP, 7 = 3(3 + n)/{b + n), is not the special one also in terms of the 
stability of the solution. As a result, the argument of the stability does not determine the 
power index of the two-point spatial correlation function in the strongly non-linear regime. 

The power index of the two-point spatial correlation function is determined only by 
the clustering process, that is, the parameter h if the self-similar solutions exist. Hence 
it is very important to estimate the parameter h and investigate whether the self-similar 
solutions exist or not in the general scale-invariant solutions which YG obtained. 
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